Zeno effect and switching of solitons in nonlinear couplers 
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The Zeno effect is investigated for soliton type pulses in a nonlinear directional coupler with dissipation. 
The effect consists in increase of the coupler transparency with increase of the dissipative losses in one of the 
arms. It is shown that localized dissipation can lead to switching of solitons between the arms. Power losses 
accompanying the switching can be fully compensated by using a combination of dissipative and active (in 
particular, parity-time -symmetric) segments. © 2011 Optical Society of America 
OCIS codes: 190.0190, 060.1810, 190.5940 



A nonlinear coupler is one of the most common and 
well studied optical devices (see e.g. 1 ). Two decades 
ago it has been realized in [2] , that a coupler whose arms 
include dissipation and gain can greatly enrich the dy- 
namical properties of the device. Remarkably, the inter- 
est in a coupler with gain and loss has recently been re- 
vitalized, since it has been linked to the so-called parity- 
time symmetric potentials [3j[4] . It has been reported 
that a passive coupler, where the dissipation acts alone 
(without gain) may also represent a significant interest. 
In particular, the phenomenon of light localization was 
observed in a linear coupler with a dissipative arm [5], 
where the phenomenon was attributed to the parity-time 
symmetry related properties. The reason for this is that 
in the linear case we can invoke a trivial transformation 
in order to reduce the pure dissipative system to a sys- 
tem with gain and loss. However, later on in [6j, it was 
argued that the increase of the transparency is intrinsic 
to the pure dissipative nature of the phenomenon when 
the dissipative losses grow. This can be observed in pure 
dissipative nonlinear systems as well. Even earlier this 
possibility was indicated in m on the basis of the anal- 
ogy between the mean-field dynamics of Bose-Einstein 
condensate and optics of Kerr media. Moreover, it was 
shown in 6 that observation of the anomalous trans- 
parency of either a linear or a nonlinear coupler has 
on the one hand relation to elementary properties of a 
mechanical pendulum with dissipation and on the other 
hand to the Zeno phenomenon. 

The Zeno phenomenon, introduced in quantum me- 
chanics in [8] and consisting in strong suppression of 
the decay of an unstable particle by means of perma- 
nent measurements, was argued [7 to have more general 
nature. This connection stems from the fact that in the 
mean-field (quasi-classical) approximation elimination of 
atoms is reduced to the effect of dissipation, which is 
described by the equations identical to those governing 
the light propagation in the unidirectional coupler [6]. 
Moreover, in 6 it was suggested how to implement a 
possible experimental setup for observation of the Zeno 
effect in optics. In particular the decay law of continuous 



wave (cw) radiation and the dependence of the damp- 
ing in the opposite arm of the coupler was analyzed as 
function of dissipation. The regimes of anomalous optical 
transparency of the coupler were suggested. These stud- 
ies, however, were concentrated on the steady beams, 
possessing no spatial-temporal structure. This raises the 
question about the possibility of observing the Zeno ef- 
fect for optical solitons, which is the main goal of the 
present Letter. 

The evolution of pulses in a nonlinear directional cou- 
pler is described by the system of coupled nonlinear 
Schrodinger (NLS) equations for the dimensionless fields 
u\ (z, t) and u 2 (z,r) in the two arms 



iUi z + Ui TT + |^l| 2 ^i + U 2 + %^\U\ = 0, 
iu 2z + U 2rr + \u 2 \ 2 U 2 + ^1 + ^72^2 = 0. 



(i) 



The variable z measures distance along the fibers and 
r measures time in a reference frame moving with the 
group velocity of the optical carrier wave. The connec- 
tion is described by coupling constants, which are set to 
be unity, as they can be incorporated into the definition 
of 7^, j — 1,2, and which are subjected to losses (or 



gain) described by 71,2 [9 10 



To outline the main idea, we first consider the propa- 
gation of solitons in each of the arms. Assuming that soli- 
ton shapes change gradually, we employ the Lagrangian 
approach in the coupler [lO], i.e., we look for solutions 
of jlj) in the form = Ajsech (r/a) e 1 ^ 3 where the 
pulse parameters, i.e. the amplitudes Ai 5 2, the duration 
a, and the phases (j>i^ are considered as slow functions 
of the propagation distance z. Using this ansatz and fol- 
lowing [To] we can calculate the averaged Lagrangian: 



L = -2A\aq. 



2Ala<j)2 



-(A( +AI) + 2aA 1 A 2 cos( 



where the amplitudes, the relative phase qb — qb\ — ob 2l and 
the pulse duration are considered as dynamical variables. 
The associated Lagrange equations lead to the following 
dynamical system (see e.g. [To] for the technical details): 
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Fig. 1. (Color online) Observation of a Zeno soliton in 
component u\ for large enough damping (72) in com- 
ponent U2 for the input pulse of the form u\{z = 0) = 
sechfr/V^), 112 (z = 0) = 0. (a) Abrupt decay of the 
field u\ for 72 = 1. (b) A Zeno soliton in u\ for 72 = 10. 
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-7(1 - F 2 ) + 2y/l - F 2 sin(0), (2a) 
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(2b) 



Qz = -71^(1 + ^-72^(1-^). (2c) 

Here F = (Pi — P2VOPL + P2) characterizes the relative 
distribution of the power P\^ (z) = |^i,2 (z,t) \ 2 dr 
between the arms and Q = (Pi + P2)/Pq is the total 
power normalized to the input power Pq at the coupler 
input Q{z = 0) = 1 and 5 = Po/3a. Notice that in the 
obtained approximation the pulse width a is constant 
which explains why Eqs. Q are identical to those ob- 
tained and studied in [6] for the coupler operating in 
the stationary regime. This also means that one can ex- 
pect various types of dissipative dynamics with solitons, 
similar to those observed in the stationary regime. 

Now we turn to the direct numerical simulations of 
Eqs. 0. First, in order to emphasize the main effects, 
we will assume that one of the arms (j = 1) is trans- 
parent, i.e., 71 = 0, and study the dynamical regimes 
depending on the losses 72 in the second arm. In Fig. [I] 
we illustrate our main result. Here we observe that ap- 
plying the input pulse to the transparent arm for a rel- 
atively weak dissipation, the output power in both arms 
display exponential decay [panel (a)], while an increase 
of the losses by a factor 10, one strongly suppresses the 
power loss in the transparent arm, giving rise to a " Zeno 
soliton" in component u\ at the output [panel (b)]. 

The two cases in Fig. [T] correspond to the input power 
below the threshold of switching of the conservative cou- 
pler (i.e. ([2| with 7^2 = 0), which for the stationary cou- 
pler is A 2 cr = 6 1 , while for the coupler operating in the 
solitonic regime it is A 2 r « 6.7 [l , 10 . 



Below the threshold the emergence of a Zeno soliton 
occurs independently of the input power Pi (0). How- 
ever, for input powers above the threshold the situation 
is changed. There appear two very different regimes of 
the soliton decay, as is shown in Fig. [2] Starting with the 
undercritical input power we clearly observe the Zeno 
phenomenon, expressed in the decrease of the decay ex- 



Fig. 2. (Color online) Distribution and decay of powers 
Pj. (a) The exponential decay in the first arm and local 
energy transfer between the arms for the initial data as 
in Fig. [l] (i.e. below the threshold, A\ < 6.7). Black solid 
curves shows the full simulation of Eq. 0, while yellow 
(gray) dashed curves shows the approximate results ob- 
tained from ([2]). (b) The same as in panel (a) but for 
the input amplitude A\ = 5, i.e. above the threshold, 
A 2 > 6.7. The horizontal dashed line shows the normal- 
ized critical threshold value A cr ~ v6?f (i.e. the value 
V / 6/7/5 c± 0.5). The sloop —2/72 is illustrated with the 
dashed-dotted line for 72 = 10. 



ponent [Fig. [2^ a)]. It is remarkably well described by the 
law Pi/Pi(0) ~ exp (—22/72). This law follows from the 
system ([2| of the adiabatic change of the soliton parame- 
ters, and hence its physics is the same as described in [6] 
for the stationary coupler [we emphasize the very accu- 
rate approximation given by the system ^ here]. The 
field amplitude in the second arm, which is zero at the 
input, achieves its maximum at short distances (being 
always a few orders below the field amplitude of the first 
component) and after that rapidly decays. We also no- 
tice that energy transfer between the transparent and 
lossy arms is strongly suppressed by the losses. If the 
input power is well above the critical value, the situa- 
tion is first inverted, and at small distances one observes 
a "standard" decay exponent proportional to 72. How- 
ever, as soon as the pulse power has decreased to some 
critical value, well estimated by A cr « 6.7 [horizontal 
dashed line in Fig. (2^b)], a change of the exponential 
decay, i.e., the Zeno effect, is observed. 

To test the accuracy of the predictions given by ([2|, 
we performed numerical comparisons of the half-width 
at the half-maximum of the solution (which for u[ s \r) 
is given by a arccosh(v // 2) used in Fig. |3[a) to normalize 
the width), and of the relative central phase <j>{r = 0, z) 
with the corresponding values of (j) obtained from ([2|, 
see Fig. [3j The peculiarity is that the width of the 
wavepacket in the lossy arm is about 1.5 times bigger 
than the width of the input soliton, which can be ex- 
plained by the dominating dispersion in the second arm 
due to low intensity of the pulse. 

It has been shown in 7 that rich possibilities of the 
control of the stationary modes in a two-mode (dimer) 
system can be achieved by using localized dissipation. 
In our case this scenario can be realized by using cou- 
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Fig. 3. (Color online) (a) The half-width at the half max- 
imum vs propagation distance for different input am- 
plitudes Ai (a = y/2/Ai) and a fixed 72 = 10. Solid 
(dashed) curves are for arm j = 1 (2). (b) The relative 
central phase <fi(z, r = 0) obtained from ([I]) (solid curves) 
and from ([2| (dashed curves), normalized with tt. 



pier arms, which are transparent everywhere except at 
localized dissipative segments. Now we turn to the dy- 
namics of solitons in such systems focusing on the possi- 
bilities of implementing different switching regimes. The 
first simple illustration of the soliton switching is shown 
in Fig. [4^a). At the input of the coupler both arms are 
transparent and the coupler is operating in the "soliton 
locked" regime (i.e. above the threshold). Without dissi- 
pation the main power would propagate in the first arm. 
If however, along the propagation distance in the first 
arm, a short dissipative segment is included, one observes 
switching of the power to the second arm (F < 0). This 
process is accompanied by the loss of more than 50% of 
the total power and with the change of the "direction" 
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of rotation of the relative phase [the inset of Fig. [4 
The switching phenomena for solitons seen in Fig. 4 
be qualitatively, but not quantitatively, reproduced by 
the approximate ODE model ([2|, which only have one 
parameter for the width a = {a\ + a^)j2. 

The problems of losses can be completely resolved 
by using simultaneously dissipative and active gain seg- 
ments. In Fig. [4^b) we use a parity-time (VT) like seg- 
ment, 72 (z) = —71(2), and observe that the total power 
do not decrease dramatically anymore. We emphasize, 
however, that in spite of the similarity of our model with 
VT structures the effect can be observed (and even opti- 
mized) using the engineering of the dissipative and gain 
defect, which do not have well defined symmetry (i.e. 
having different lengths or amplitudes). Finally, in the 
inset of Fig. [4^b), we have engineered the amplitudes 
(r:s) and the translations (z a y.s) of the first dissipa- 
tion/gain segment, and in addition added a second seg- 
ment of gain/dissipation in order to switch back, and 
with the total power being kept close to Q(z) « 1. 

To conclude we have reported the possibility of obser- 
vation of the Zeno effect with solitons in a coupler with 
one transparent and one lossy arm. The effect consists in 
the increase of the total transparency of the coupler when 
the losses of one of the arms are increased. We have also 
shown that using dissipative segments one can perform 
efficient management of solitons, in particular switching 
them at will. Since the intensity losses are inevitable in 




Fig. 4. (Color online) Switching between the arms in- 
duced by localized defects. In (a) only a dissipative seg- 
ment 71(2) = r (arctan[5(z — z a )] — arctan[5(z — 2^)] ) 2 
with T = 0.135, z a = 1.5, and z^ = 3.0 in the first arm is 
included. The inset shows the relative phase. The solid 
(blue) and dashed (red) curves show the simulations of 
the models in ([!]) and respectively. In (b) the simul- 
taneous effect of the dissipative and active segments, rep- 
resenting a VT-like structure, i.e. 72(2) = — 7i(z), with 
71(2) as in panel (a) but with T = 0.065. For the inset 
of (b) we have added a segment of gain/dissipation at a 
Az ~ 10 distance, and used engineered parameters in or- 
der to keep Q(z) « 1. In both (a) and (b) we use u\(z = 
0) = 20sech(10 v / 2r) and u 2 (z = 0) = 5 sech(5r/ v / 2), 
i.e. initially F(0) = 3/5. 



the described systems, we have suggested simultaneous 
use of the dissipative and gain segments, in order to per- 
form lossless manipulations of the solitons. 
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